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ABSTRACT 

We here present two-dimensional, time-dependent radiatively cooling hydrodynamical simu- 
lations of the large and little Homunculus nebulae around rj Carinae. We employ an alternative 
scenario to previous interacting stellar wind models which is supported by both theoretical and 
observational evidence, where a non-spherical outburst wind (with a latitudinal velocity depen- 
dence that matches the observations of the large Homunculus), which is expelled for 20 years, 
interacts with a pre-eruptive slow wind also with a toroidal density distribution, but with a 
much smaller equator-to-polar density contrast than that assumed in previous models. A sec- 
ond eruptive wind with spherical shape is ejected about 50 years after the first outburst, and 
causes the development of the little internal nebula. We find that, as a result of an appropriate 
combination of the parameters that control the degree of asymmetry of the interacting winds, 
we are able to produce not only the structure and kinematics of both Homunculus, but also the 
high-velocity equatorial ejecta. These arise from the impact between the non-spherical outburst 
and the pre-outburst winds in the equatorial plane. 
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The dusty large (16" long) bipolar Homunculus 
ejecta around rj Car is a hollow reflection nebula 
(e.g., Smith et al. 2003a and references therein) 
that was produced by the 20 year Great Eruption 
of the star, from ~1840 to 1860. Its axis is in- 
clined ~ 45° to the line of sight (Davidson et al. 
2001), and at a distance of ~2.3 kpc, it has a total 
physical size ~ 6 x 10 17 cm. The hot features ob- 
served outside the bipolar nebula were very prob- 
ably ejected earlier in episodic events before the 
major eruption (Walborn et al. 1978; Weis et al. 
2001). Recently, an inner bipolar emission neb- 
ula has been also discovered embedded within the 
larger Homunculus (extending from —2" to +2" 
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across the star) which may have been originated 
from a minor eruption event in the 1890s, i.e., ^50 
years after the formation of the larger Homuncu- 
lus (Ishibashi et al. 2003). This little Homunculus 
seems to follow approximately the shape of the 
larger one. Observations also evidence the exis- 
tence of ejecta in the equatorial region that may 
contain material from both the 1890 eruption and 
the great eruption in the 1840s (Davidson et al. 
2001). The velocity of this material obtained from 
Ha profiles may reach velocities ~ 400 — 750 km 
s" 1 (Smith et al. 2003a). 

As an extreme luminous blue variable star 
(LBV), -q Car loses copious amounts of mass in 
form of quasi-steady winds punctuated by erup- 
tive events where the mass loss may increase by 
at least an order of magnitude in short periods of 
time (Maeder 1989; Pasquali et al. 1997). Cur- 
rently, the dominant velocities in the expanding 
Homunculus are 400 to 600 km s _1 , but some 
faster material (~ 1000 km s _1 ) has been de- 
tected in the poles (e.g., Smith et al. 2003a), 
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and the mass-loss rate is ~ 10 4 — 10 3 M Q yr 1 
(Humphreys & Davidson 1994, Hillier et al. 2001; 
Corcoran et al. 2001; Soker 2001). The inner 
bipolar nebula, produced during the 1890s erup- 
tion, has a peak velocity ~ 300 km s _1 and a total 
estimated mass ~ 0.1 M Q (Ishibashi et al. 2003). 

It has been previously suggested that the shap- 
ing of the 77 Car nebulae could be explained by a 
colliding wind binary star model. Soker (2001), 
for example, has argued that the companion star 
could divert the wind blown by the primary star, 
by accreting from the wind and by blowing its 
own collimated fast wind that could have, in turn, 
played a role in the formation of the Homuncu- 
lus lobes. A strong argument against a compan- 
ion star dominating the wind structure is that it 
appears to be predominantly symmetric, as indi- 
cated by recent STIS spectral observations in sev- 
eral positions along the Homunculus (Smith et al. 
2003a). These show the same latitudinal depen- 
dence for the velocities in both hemispheres and 
both sides of the polar axis, and the same P Cygni 
absorption in hydrogen lines on either side of the 
poles. Although a colliding wind binary model 
cannot be disregarded at the present, we here as- 
sume that the shaping of 77 Car nebulae is domi- 
nated by the primary star's wind. 

As noticed by Smith et al. (2003a), the high ve- 
locities seen in reflected light from the polar lobes 
give a first direct evidence that the polar axis of 
the Homunculus is aligned with the rotation axis 
of the central star. This has important conse- 
quences for the formation of the bipolar lobes and 
the equatorial ejecta around 77 Car, as it may be 
an indication that axial symmetry and the ejec- 
tion mechanism during the Great Eruption were 
directly linked to the central star's rotation. Also, 
the observed latitudinal variations in H and Hcl 
lines revealing that the speed, density and ioniza- 
tion in 77 Car wind are non-spherical nearby the 
star, may be an indication that the stellar wind is 
inherently non-spherical. 

In previous work, Frank et al. (1995) have 
performed low resolution two-dimensional numer- 
ical simulations of the large Homunculus of 77 
Car, adopting an interacting stellar wind scenario 
wherein a spherical fast wind expands into a non- 
spherical (toroidal) slow, dense wind previously 
ejected from the star. They found that the Ho- 
munculus morphology could be reproduced with 



an equator-to-polar density ratio ~ 200, which 
would imply the existence of a very dense toroidal 
environment surrounding the nebula. Langer et 
al. (1999), have assumed a variant of this scenario 
including the effects of stellar rotation. Using the 
wind-compressed model of Bjorkman & Cassinelli 
(1992), they showed that a strong equator-to-pole 
density contrast could have formed during the 
great outburst in the 1840s if the star was close 
to the Eddington luminosity limit. Under this 
circumstance, the centrifugal and radiative forces 
must balance gravity at the equator, and a strong 
non-spherical mass loss should occur deflecting the 
wind streamlines towards the equator. In this 
model, a spherical, fast post-outburst wind pro- 
duces the bipolar bubble through interaction with 
the toroidal outburst flow. In addition to form- 
ing lobes with an approximate shape to the ob- 
served Homunculus, this numerical model (which 
has included the effects of time-dependent radia- 
tive cooling) has revealed the development of small 
fingers at the shell surface caused by Vishniac type 
instabilities. 

Although these previous models are partially 
successful at reproducing the basic shape of the 
large Homunculus nebula, they both rely on the 
presence of a thick torus around the Homuncu- 
lus nebula (which is very dense in the equato- 
rial region). Observations however, indicate only 
the presence of a faint nebulosity surrounding the 
large Homunculus. To overcome this difficulty, 
more recently Frank et al. (1998), and Dwarkadas 
& Balick (1999) have proposed alternative mod- 
els. Assuming an inverted scenario, in which a 
non — spherical fast wind expands into a previ- 
ously deposited isotropic slow wind, Frank et al. 
(1998) have found that they are able to reproduce 
strongly bipolar outflows with polar caps that can 
be denser than the lobes' flanks. However, like the 
previous ones, this model is unable to produce the 
equatorial ejecta. Dwarkadas & Balick (1998), on 
the other hand, have replaced the thick torus of 
the previous models by a small and dense, near- 
nuclear toroidal ring. This also manages to pro- 
vide some collimation of the spherical wind ejected 
during the Great Eruption. Besides, in the pres- 
ence of radiative cooling, the ring is completely 
destroyed by the impact of the wind, and the au- 
thors have claimed that this fragmentation of the 
ring could help to explain the equatorial ejecta. A 
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potential problem with this interpretation is that 
their model predicts velocities for the fragments 
that are too small (~ 50 — 100 km s _1 ) compared 
to the observed ones in the outer parts of the ejecta 
(e.g., Smith et al. 2003a). 

A successful modeling of the formation of the 
large and little Homunculus nebulae around r\ Car 
should account for both the bipolar morphology 
and the equatorial ejecta. We here present results 
of numerical simulations that consider an alterna- 
tive scenario to the interacting stellar winds mod- 
els above, in which a fast, non-spherical wind is 
ejected for 20 years (with a latitudinal velocity 
dependence that matches the observations in the 
large Homunculus), interacts with a pre-outburst 
slow wind also with a toroidal density distribution, 
but with a much smaller equator-to-polar density 
contrast than that assumed in previous models. 
A second eruptive wind with spherical shape is 
ejected about 50 years after the first outburst, and 
causes the development of the little internal neb- 
ula. We find that, as a result of an appropriate 
combination of the parameters that control the 
degree of asymmetry of the interacting winds, this 
model is able to produce not only the structure of 
both Homunculus nebulae, but also the equatorial 
ejecta (see below). 3 

We have carried out several hydrodynamical 
simulations considering various possible scenarios 
for the degree of asymmetry of the interacting 
winds, but here, we will present only the model 
that has best matched the observations. A more 
detailed description of the other models will be 
presented in a forthcoming paper (Gonzalez, de 
Gouveia Dal Pino, Raga & Velazquez 2003). In 
contrast to the previous works, our simulations 
compute explicitly the time dependent radiative 
cooling of the gas including several atomic and 
ionic species, which allow for a more realistic eval- 
uation of its effects on the flow. 



We note that there is a number of proposed mechanisms in 
the literature that predict the development of intrinsically 
non-spherical winds coming out from rotating LBV stars 
(see, e.g., Bjorkman & Cassinclli 1992; Lamers & Pauldrach 
1991; Owocki et al. 1996, 1998). 



2. Numerical Method, Initial Physical 
Conditions, and Results 

In order to investigate the winds' interaction, 
we have performed gasdynamic simulations using 
a modified two-dimensional version of the Yguazu- 
a adaptive grid code originally developed by Raga 
et al. (2000; see also Raga et al. 2002, Masciadri 
ctal. 2002, Velazquez et al. 2003). This code inte- 
grates the hydrodynamic equations explicitly ac- 
counting for the radiative cooling together with a 
set of continuity equations for several atomic/ionic 
species employing the flux-vector splitting algo- 
rithm of Van Leer (1982). The following species 
have been considered: H I, H II, He I, He II, He 
III, C II, C III, C IV, O I, O II, and O III. The 
calculations were performed on a five-level, binary 
adaptive grid with a maximum resolution along 
the x and y axes of 7.81 x 10 14 cm. The computa- 
tional domain extends over (4 x 10 17 ) x (4 x 10 17 ) 
cm, corresponding to 512 x 512 grid points at the 
highest resolution grid level. 

We assume that a light, hot gaseous toroidal 
distribution was formed around i] Car prior to the 
Great Eruption of -1840 (e.g., Weis 2001). To 
produce it, a slow and steady wind is assumed to 
emanate from the stellar surface into the ambient 
medium (with initial temperature T a = 100 K and 
number density n a = 10~ 3 cm~ 3 ), with the follow- 
ing properties (see Frank et al. 1995), 

v = voF , (2) 

where tq is the injection radius, vq and no are 
the velocity and the number density at the pole, 
respectively. We have considered a similar func- 
tion to that adopted by Frank et al. (1995, 1998) 
to produce a smooth variation in the wind den- 
sity and velocity from the equator to the pole, 
Fg = 1 - a[(l - e -2^n 2 e )/(1 _ e -2/3^ where q 

is the polar angle, the parameter (3 controls the 
shape of the wind and a the pole-to-equator den- 
sity and velocity contrasts. The density and ve- 
locity angular dependence in equations (1) and 
(2), respectively, imply a constant mass- loss rate 
(M tx n v) as a function of the angle. 

We have assumed for the pre-eruptive wind a 
mass loss rate M = 10~ 3 M Q yr" 1 ; v = 250 km 
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s -1 ; a — 0.9; (3 — 1.5; r = 1 x 10 16 cm; and ini- 
tial temperature To = 10 4 K. When the slow wind 
reaches the edge of the computational domain in 
the polar direction (y-axis), a non-spherical out- 
burst wind is turned on, for 20 years, to produce 
the Great Eruption, with M ~ 7 x 10~ 2 M Q yr _1 ; 
vo = 715 km s _1 ; and To = 10 4 K. For simplic- 
ity, we have adopted the same functional depen- 
dence [F(0)] as above to compute the velocity and 
density angular variations in this wind, but with 
different values for a and (3 (0.78 and 0.3, respec- 
tively). These values were obtained from the best 
fit of equation (2) to the observed latitudinal vari- 
ation of the expansion velocity of the large Ho- 
munculus (Smith 2002; Davidson ct al. 2001). 
After 20 years, a third wind with the same con- 
ditions of the original slow wind resumes for 30 
years and then, another outburst (which was as- 
sumed, for simplicity, to be spherical) is allowed 
to occur for about 10 years with M = 10~ 2 M 
yr -1 ; and v — 317 km s _1 , after which the origi- 
nal slow wind again resumes. The adopted values 
of a and (3 for the slow pre-outburst wind imply a 
much smaller (larger) equator-to-pole density (ve- 
locity) contrast (= 10) than that used by Frank 
ct al. (1995). This results a fainter and lighter 
toroidal envelope around the large Homunculus, 
as required by the observations. 

Figures 1 and 2 depict the results of this simula- 
tion with the interaction of the five winds above. 
We notice that two bipolar expanding shells de- 
velop both with shapes and kinematics similar to 
the large and the little Homunculus. As in pre- 
vious calculations (e.g., Frank et al. 1995, 1998), 
we find that during the 160 years of evolution, the 
momentum flux is dominated by the outbursts, so 
that the post-outburst slow winds have no signif- 
icant effect on the formation of the Homunculus 
structures. In the case of the outer Homunculus, 
since its mass-loss rate is almost an order of mag- 
nitude larger than that of the pre-outburst wind, it 
initially expands almost ballistically (in the polar 
direction), without being much decelerated by the 
pre-outburst slow wind, and retains most of the 
non-sphericity imprinted in it near the star, as ob- 
served. In the equatorial direction, this first out- 
burst wind reaches the pre-eruptive toroid when 
its shock front is still within the computational 
domain. The impact of the two fronts (at t ~ 100 
yr after the great eruption; see Fig. lc) causes the 



formation of a faint equatorial ejection that could 
explain the outer parts of the observed equatorial 
skirt in r\ Car (Fig. Id). This ejection moves with 
a mean velocity ~ 700 km s _1 , therefore in quali- 
tative agreement with the observations (Smith et 
al. 2003a). 

The internal bubble that develops from the im- 
pact of the second (spherical) outburst with a pre- 
outburst wind ejected with the same characteris- 
tics of the first slow wind (Fig. lb), soon acquires 
a bipolar shape (Fig. lc) which is very similar 
to that implied from the observations of the little 
Homunculus with an expansion velocity ~ 300 km 
s" 1 (Ishibashi et al. 2003). 

The analysis of the results at the outer Ho- 
munculus shows that a double-shock structure has 
developed with an outward moving shock that 
sweeps the material of the precursor wind and an 
inward shock that decelerates the outburst wind 
material coming behind. A thin dense and cold 
shell develops at the edge of the nebula as a result 
of the strong compression that follows the radia- 
tive cooling of the shocked material behind both 
shocks. For a strong shock (with v s > 80 km s _1 ), 
the radiative cooling time for an one-dimensional 
shock is given by t c ~ 320yr u^iqo Ppre ,-22' where 
v s ,ioo is the shock speed in units of 100 km s _1 
and pp re is the pre-shock densitiy in units of 10~ 22 
g cm -3 (Hartigan, Raymond & Hartmann 1987; 
Gonzalez 2002). At the poles, the inward shock 
speed in the simulation is v c ~ 80 km s _1 , and 
p pre ~ 5.1 x 10~ 20 g cm~ 3 . This implies t c ~ 0.5 
yr, which is very short compared to the age of the 
large Homunculus bubble (~ 160 yr), and results 
in very small cooling distance behind the inner 
shock that can be estimated from d c ~ 3.7 x 10 14 
cm w 4 ioo Ppre, -22 (Hartigan, Raymond & Hart- 
mann 1987), or d c ~ 2 x 10 11 cm behind the in- 
ward shock. This value qualitatively explains the 
narrowness of the cold thin shell seen in the sim- 
ulations (Fig. 1), which, at t — 160 yr, has a tem- 
perature ~ 6000 K and a density 3 x 10~ 19 g cm -3 . 
Behind the outward shock, on the other hand, the 
smaller pre-shock density (p pre - 2.1 x 10~ 21 g 
cm~ 3 ) and the higher shock speed (v c ~ 400 km 
s _1 ) produce a larger cooling time (t c ~ 69 yr) 
and, as a result, a thicker and hotter polar cap 
with d c ~ 10 16 cm (according to the equation 
above), and T ~ 1.6 x 10 4 K and p ~ 2.1 x 10~ 21 g 
cm~ 3 (as obtained from the simulations in Fig. 
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1). Despite the approximations involved in the 
evaluation of the cooling distance above, it is only 
a factor three smaller than the thickness of the 
outer shell obtained from the simulations which is 
<~ 3.5 x 10 16 cm. This value is in turn, comparable 
to the observations that indicate a radial thickness 
of the outer Homunculus ~ 1" = 2500 AU at the 
poles (Smith et al. 2003a). 

A similar analysis for the expanding shell of the 
inner Homunculus shows that its thickness at the 
poles (obtained from the simulations) is ~10 16 cm. 

3. Discussion and Conclusions 

The results of our 2-D hydrodynamical simula- 
tions involving the interaction of five winds with 
different initial conditions indicate that the shape 
and kinematics of the large Homunculus of n Car 
can result from the interaction between fast and 
slow intrinsically non — spherical winds. This 
model is a variant of previous interacting wind 
scenarios that have assumed either fast spherical 
winds interacting with a dense and heavy toroidal 
environment (Frank et al. 1995; Langer et al. 
1999; Dwarkadas & Balick 1999) or fast non- 
spherical winds interacting with an isotropic envi- 
ronment (Frank et al. 1998). It has two attractive 
advantages: (i) it shows that a non-spherical fast 
wind impinging on a slow toroidal wind is able to 
produce the high- velocity outer parts of the equa- 
torial ejecta observed around r\ Car (Smith et al. 
2003a) ; and (ii) the choice of a lighter pre-outburst 
wind in our model, has resulted in a less dense 
toroidal halo around the large Homunculus neb- 
ula than in previous models, as required by the 
observations. 4 

It is noteworthy that, in numerical experi- 
ments where, instead of an initially non-spherical, 
a spherical outburst wind was injected into a 
toroidal pre-outburst slow wind with the asym- 
metry parameters a and (3 determined either from 

4 We note that the equatorial waist produced in the sim- 
ulations is slightly thicker than the one observed. While 
observations indicate that the ejections are <0.4 times the 
diameter of the Homunculus lobes, in the simulations this 
is ~0.6. We also note that the ratio of the length to the 
width of the two Homunculus is somewhat larger in the 
simulations. Both differences could be partially attributed 
to projection effects which were not considered here, and 
could be diminished with minor changes in the adopted 
parameters for the model. 



the observed large Homunculus expansion veloc- 
ity distribution or given by the same values as 
those used by Frank et al. (1995) for a heavier 
toroid, have failed to produce simultaneously the 
shape of the large Homunculus and the equato- 
rial ejecta. In these cases, the encounter of the 
shock fronts of the two winds first in the equato- 
rial region and afterwards in the lateral regions 
of the outer bubble causes its fragmentation and 
spreading of the shell material at high latitudes, 
thus destroying the bipolar morphology (see Gon- 
zalez et al. 2003). These results suggest that, in 
order to produce both the Homunculus bipolar 
morphology and the equatorial ejection from the 
winds interaction, these must be both intrinsically 
non-spherical (as in Figs. 1 and 2). 

Although the interaction of the second out- 
burst wind (assumed to be spherical) with its pre- 
outburst wind was able to produce the internal 
Homunculus, it has failed to develop internal equa- 
torial ejecta in the simulation depicted in Figures 
1 and 2. However, an appropriate combination of 
non-spherical wind parameters in this case simi- 
lar to the one of the 1840s outburst, could also 
probably generate an internal equatorial ejection. 
In fact, recent UV images within 0.2 arcsec of the 
star, have revealed the existence of a little internal 
torus that may be related to the little Homuncu- 
lus and may signify that a recurrent mass ejection 
with the same geometry as that of the Great Erup- 
tion may have occurred (Smith et al. 2003b, see 
also Gonzalez et al. 2003). 

In order to simulate an experiment using a con- 
dition at the base of the r/ Car wind similar to 
the one presently suggested by the observations 
(Smith et al. 2003a), we have also computed a 
model in which the non-spherical outburst wind of 
1840s impinges on a slow pre-outburst wind with 
a larger density (and mass-loss rate) in the polar 
direction (n oc Fg in cq. [1]). We find that this 
scenario is unable to develop a narrow equatorial 
ejecta. This result suggests that the conditions at 
the wind base prior to the Great Eruption in the 
1840s were probably not the same as the current 
ones. 

Finally, we notice that, despite the high- 
resolution and the explicit time-dependent compu- 
tation of the radiative cooling of the gas, our sim- 
ulations, similarly to the radiative cooling models 
of Frank et al. (1998), have not revealed the for- 
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mation of small fragments on the surface of the 
Homunculus, as those seen in Langer et al. (1999) 
simulations, which have resulted from Vishniac in- 
stabilities in the radiativcly cooled shell long after 
the eruption. This is probably due to differences 
in the initial conditions between the two models. 
Nonetheless, a granular structure is effectively ob- 
served on the large Homunculus surface. Is is not 
improbable, however, that they have resulted from 
variability or instabilities in the winds near the 
surface of the star (Smith ct al. 2003a). This 
question, as well as three-dimensional effects, will 
be addressed in future work. 
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Fig. 1. — Gray-scale map of the density distribution 
(in logio scale) for four different times in the evolution 
of the model, (a): t = 10 yr; (b) t = 60 yr; (c) t = 110 
yr; and (d) t = 160 yr. The density (at the vertical 
scale on the right side of the figures) is in g cm~ 3 , and 
the x and y axis are in cm. 



Fig. 2. — The same as in Figure Id where it is de- 
picted the present-day structure of the system with the 
inner and outer Homunculus nebulae, and the equato- 
rial ejections. The density map has been rotated by 
an angle of 45°, like in the observations. 
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